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Abstract 

We study five dimensional thin-shell wormholes in Einstein-Maxwell theory with a Gauss- 
Bonnet term. The linearized stability under radial perturbations and the amount of exotic 
matter are analyzed as a function of the parameters of the model. We find that the inclusion 
of the quadratic correction substantially widens the range of possible stable configurations, and 
besides it allows for a reduction of the exotic matter required to construct the wormholes. 
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1 Introduction 

Traversable Lorentzian wormholes are solutions of the equations of gravitation which connect two 
regions of the same universe, or of two universes, by a throat m i2]- The throat is defined as a 
minimal area hypersurface which satisfies a flare-out condition [51- This requires the presence of 
exotic matter, that is, matter which violates the null energy condition (NEC) jT] 2JQL 3j. Because 
it has been shown that the amount of exotic matter necessary for the existence of a wormhole can 
be made infinitesimally small by a suitable choice of the geometry [5j, considerable efforts have 
been addressed to precisely quantify such amount for different configurations, and to show how it 
can be minimized j||. Indeed, the total amount of exotic matter has been pointed as an indicator 
of the physical viability of traversable wormholes jjj. 

Another central aspect of a wormhole -in fact, of any physically meaningful solution within any 
theory of gravitation- is its stability. Within the framework of traversable wormholes, stability un¬ 
der perturbations preserving the symmetry of the original configuration has been widely analyzed. 
In particular, this problem has received considerable attention in the case of thin-shell wormholes, 
that is, wormholes which are mathematically constructed by cutting and pasting two manifolds to 
obtain a geodesically complete new manifold 0. In this case the exotic matter is located in a shell 
placed at the joining surface; the framework for dealing with these wormholes is the Darmois-Israel 
formalism, which leads to the Lanczos equations, that is, to the equations of gravitation projected 
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on the joining surface HUD]. The solution of these equations, once provided an equation of state 
for the matter on the shell, determines the dynamical evolution. Such a procedure has been applied 
to spherically and cylindrically symmetric configurations associated to wormhole solutions within 
general relativity lT. I2L fl~3l. 

The theory of gravity in five dimensions corresponding to the Einstein-Hilbert action supple¬ 
mented with a Gauss-Bonnet term is, in a certain sense, the most general (metric) theory of gravity 
one can construct satisfying the conservation of the equations of motion which still remain being 
of second order M- This theory, and its analogue in D dimensions, was extensively studied in 
the last three decades and, in particular, the attention focused on it was mainly due to the fact 
that the theory arises within the string theoretical framework m- For instance, a version of this 
appears as corrections (proportional to the inverse of the string tension) to the low energy effective 
action of the heterotic string theory m as well as in Calabi-Yau compactification of the M-theory 
m- Besides, these theories result closely related to Chern-Simons gravity in odd dimensions which 
turns out to correspond to a particular choice of the parameters of the model; see Ref. m Physi¬ 
cally, adding the Gauss-Bonnet (higher order) terms in the gravitational action corresponds to the 
inclusion of short distance corrections to general relativity. The study of black hole solutions in 
Einstein-Gauss-Bonnet theory was initiated in the decade of 1980, when the statical spherically 
symmetric solution was reported by Boulware and Deser in Ref. m- Subsequently, Wiltshire 
derived in Ref. m the charged black hole geometry in both Maxwell and Born-Infeld electrody¬ 
namics. Both geometrical and thermodynamical aspects of black holes are substantially modified 
by the addition of the Gauss-Bonnet term m- 

Lorentzian wormholes in spacetimes with more than four dimensions were analyzed by several 
researchers Esina. In particular, wormholes in Einstein-Gauss-Bonnet gravity were considered 
in Ref. 122 . In the present work, the idea is then to use Lorentzian wormholes as a test bed to 
explore some of the qualitative changes that could happen in General Relativity with the addition 
of a Gauss-Bonnet term. More precisely, we shall study how the stability under radial perturbations 
of five dimensional spherically symmetric thin-shell wormholes in Einstein-Maxwell theory and the 
amount of exotic matter needed are affected by the presence of a Gauss-Bonnet term. In Section 
2 we shall construct a generic thin-shell wormhole and write down the corresponding Lanczos 
equations; in Section 3 we shall analyze its mechanical stability under perturbations preserving the 
symmetry, and in Section 4 the energy conditions will be studied, and the total amount of exotic 
matter will be calculated. The dependence of the results in terms of the parameters of the model 
will be analyzed in detail. It will be shown that the inclusion of the Gauss-Bonnet term permits 
stability configurations with more physical values of /3 2 with small charge; also, we will see that the 
amount of exotic matter can be reduced for given values of the parameters. Section 5 is devoted 
to a brief summary and discussion. Throughout the paper we set units so that c = G = 1. 


2 Charged thin-shell wormholes 


The five-dimensional Einstein-Maxwell theory with a Gauss-Bonnet term representing a quadratic 
curvature correction is given by QHH] 


s 


d 5 Xy/—g 


R- 2A 


\f^F^ + a ( R a p lS R a ‘^ - 4 R a0 R a P + R 2 ) , 


( 1 ) 


where the signature chosen is (—b + + +), A is the cosmological constant, and a is a constant of 
dimensions (length) 2 . The variational principle 6S = 0 leads to the Einstein-Maxwell equations 
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where is the usual electromagnetic energy-momentum tensor, and TpP is an effective tensor 
associated with the quadratic Gauss-Bonnet term included in the action. These equations admit 
a spherically symmetric solution given by HI 
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The non null components of the electromagnetic tensor in an orthonormal frame are F;~ = —F-a = 


Q/Anr 3 . It is not difficult to check that in the limit a —> 0 the five dimensional Einstein-Maxwell 
solution with cosmological constant is recovered. In this limit, for A = 0 the five dimensional 
Reissner-Nordstrom metric is obtained, so M > 0 and Q can be identified with the mass and 
charge respectively. For a p 0, there is a minimum value of the radial coordinate r m .; n such 
that the function inside the square root in Eq. © is positive for r > r m i n , so the metric © is 
well defined. The geometry has a curvature singularity at the surface defined by r = r min M- 
Depending on the values of the parameters, this singular surface can be surrounded by an event 
horizon with radius so the metric © represents a black hole, or in absence of the event horizon, 
it is a naked singularity. 

We shall construct a spherically symmetric thin-shell wormhole starting from the generic ge¬ 
ometry ©, and introduce the explicit form © in the final results. We take two copies of the 
region r > a, with a greater than r^ and r m i n to avoid possible horizons and singularities in our 
geometry, and paste them to obtain a geodesically complete new manifold with a matter shell at 
the surface r = a, where the throat of the wormhole is located. The procedure follows the steps of 
the Darmois-Israel formalism; in terms of the original coordinates X 7 = (t, r, 9, x, ip), on the shell 
we define the coordinates £* = (r,0, %,</?), with r the proper time. Thus, using an orthonormal 
basis {ef, es, e%, e^}, the extrinsic curvature at the two sides of the shell reads 
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where nP are the unit normals to the surface. Defining k~- = K -1 — K~~ and k = tr(n~~), we obtain 

i ij ij v 

the Lanczos equations (the Einstein’s equations on the shell) 


K ij - K 9~ij = - §7r % ( 8 ) 

where = diag(—1,1,1,1) and SA = diag (cr,p§,Px^Pp is the surface energy-momentum tensor. 
To allow for the analysis of radial perturbations, we let the throat radius to vary with the proper 
time: a = a(r). As a consequence of the generalized Birkhoff theorem proved in Ref. EDI, the 
geometry will remain given by © and © for any r greater than a(r). The resulting expressions 
for the energy density and pressures for a generic metric function / turn to be 

S Tr = a = + d 2 , (9) 
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where the overdot and the prime means, respectively, the derivatives with respect to r and r. As it 
was to be expected, the energy density is negative, revealing the existence of exotic matter at the 
shell. It is easy to see from Eqs. © and m that the following conservation equation is fulfilled: 
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For a static configuration of radius ao we simply have 
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Note that as the wormhole radius approaches which would be the event horizon radius in the 
original metric, the energy density approaches to zero, but, instead, the pressure diverges unless 
f'( r h ) = 0; this will be discussed in detail in section 4. 


3 Stability analysis 

We shall study the stability of the configuration under small perturbations preserving the symmetry; 
for this we shall follow the procedure first applied to thin-shell wormholes in Ref. EU The 
dynamics of the shell results from the Eqs. © and m or, alternatively, from one of them and 
the conservation equation (HU)- In any case an equation of state for the matter on the shell must 
be provided; because we are interested in studying small radial perturbations around a radius of 
equilibrium ao, we propose a linear relation 

p = p 0 + /3 2 (cr - a 0 ), (13) 


where ao and po are the energy density and pressure corresponding to the equilibrium radius of the 
wormhole; for ordinary matter f3 2 would be the speed of sound, but due to the presence of exotic 
matter we shall regard it as an arbitrary constant. 

These equations lead to an explicit relation between the energy density and the radius: 
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Introducing this in Eq. © we obtain the equation of motion 


a 2 + V{a) = 0, 


(15) 


where the potential V(a) is defined as 

V( a ) = f{a)-^-a 2 cr 2 . (16) 

It is easy to verify that the potential fulfils V(ao) = V'(ao) = 0, so that the stable equilibrium 
configurations correspond to the condition V”(ao) > 0. The second derivative of the potential for 
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the generic metric © is given by 
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Hence the stability of the configuration requires the relation 
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if /'(ao) — 2/(ao)/ao < 0. Although it is possible to write explicitly an analytic expression of f3 2 
as a function of the parameters, the complexity of the formulas m and inhibits to get a 
clear insight of the stability regions. It is thus preferable, as it is customary, to draw the curves 
V"(ao) = 0, which allow for an intuitive understanding of the behavior of the configuration. We 
shall consider the cases A = 0 and A / 0 separately, and thoroughly discuss the dependence of the 
regions of stability with different choices of the parameters. 


3.1 Case A = 0 


As pointed above, the square root of the parameter a. of the theory introduces a length scale 
such that the Gauss-Bonnet corrections become relevant when the typical dimensions of a given 
configuration are of order ^/|a|. Thus, it is natural to perform an analysis for different cases 
depending on the relation existing between a and the mass M. 

In the case of null cosmological constant the metric from which we start presents a singularity 
located at r m i n given by the greatest real and positive solution of the equation 

6 16 Ma 2 8 Q 2 a 

r H- r -= 0. 
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If Eq. (1201) has no real positive solutions we have r m i n = 0, where the metric diverges. This 
singularity is surrounded, in principle, by two horizons with radii 
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The event horizon is placed at = r+, and r_ is the inner horizon. For a > —M/ir, r m i n < rh 
and the singularity can be shielded by the event horizon. But when a < —M/n, we have a naked 
singularity because r m i n > r^. For \a\ < M/n, it is easy to see that there exists a critical value of 
the charge 

\Q C \ = V3 


M 


— a 


7r 


( 22 ) 


5 




















Figure 1: The stability regions (marked with an S) are shown for A = 0 and a/M = 0.1, which 
implies \Q C \/M = 0.38. Only physically admissible regions, r > r/, or r > r m i„, are considered (see 
the text). For the upper left figure \Q\ = 0; in the upper right figure, \Q\ = 0.99|Q C |; in the lower 
left, \Q\ = 1.01|Q C |; and for the last one \Q\ = 2\Q C \. 


such that if \Q\ < \Q C \ there are two horizons, if \Q\ = \Q C \ there is only one (degenerate) horizon, 
and if | Q | > \Q C \ there are no horizons. For |a| > M/n, no horizons exist for any value of the 
charge (except the nonphysical solution = 0 for a = M/tt and Q = 0). As mentioned in Section 
2, the wormhole radius ao is taken greater than r/, to avoid the presence of event horizons in our 
wormhole geometries. Note that when there is no horizon, the presence of the singular surface in 
r = r m i n compels us to consider values of ao greater than r m * n for the radius of the wormhole 
throat. 

The associated stability analysis reveals both analogies and remarkable differences with the 
general relativity case corresponding to a —> 0: 1) Two distinct regimes take place for |a| < M/n 
and |a| > M/-k. In the first range most relevant results regarding the stability of the solutions 
appear, while in the second one stability requires /3 2 < 0 (Figs. HJandEJ). 2) For a > 0 the critical 
value of charge, \Q C \, is smaller than the value corresponding to general relativity; instead for 
a < 0 the critical charge is larger than in the absence of the Gauss-Bonnet terms. For a —> 0 (five 
dimensional Reissner-Nordstrom metric) the critical value of the charge is \Q C \ = 3) As 

larger values of a are considered, the regions of stability become enlarged (including (3 2 < 1, which 
is an interesting feature) without the necessity of large values of the charge; in particular (see upper 
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Figure 2: The stability regions are shown for A = 0 and a/M = 0.25, which implies \Q C \/M = 0.12. 
For the upper left figure \Q\ = 0; in the upper right figure, \Q\ = 0.99|Q C |; in the lower left, 
\Q\ = 1.01|Q C |; and for the last one \Q\ = 2\Q C \. 


left in Fig. 2), this is possible with zero charge, which constitutes a drastic difference comparing 
with the general relativity case. Besides, for \Q\ ^ \Q C \ we find a range of radii ao for which stability 
is achieved with any value of the parameter f3 2 (lower left in Fig. 2). 4) For —M/tt < a < 0 the 
regions of stability turn to be smaller than without the Gauss-Bonnet quadratic contribution (see 
Fig. 3). 


3.2 Anti-De Sitter case 

The presence of the cosmological constant introduces a restriction on the admissible values for or. 
it must be Aa > —3/4, in order to keep the metric real for large values of r. Now the singular 
surface radius r m i n is given by the greatest real and positive solution of the equation 

4AaA 16 Mol 2 8 Q 2 a . . 

1 + -5“ ) r +- r - — = °> ( 23 ) 

o J 7T O 

or r m i n = 0 if all real solutions are non positive. The addition of the cosmological constant makes 
more complicated the structure of the horizons in the original manifold. Indeed we now have that 
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Figure 3: The stability regions are shown for A = 0 and a/M = —0.1, which implies \Q C \/M = 
0.725. For the upper left figure \Q\ = 0; in the upper right figure, \Q\ = 0.99|Q C |; in the lower left, 
\Q\ = 1.01|Q C |; and for the last one \Q\ = 2\Q C \. 
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Vm Vm 


Figure 4: The stability regions are shown for A = 0 and a/M = 1/ir, which implies \Q C \/M = 0. 
The left figure corresponds to \Q\/M = 0, and the right one to \Q\/M = 0.1. 


the horizons are real and non-negative solutions of the equation 
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Figure 5: The stability regions are shown for A = 0 and a/M = 1.01 /7r. The left figure corresponds 
to \Q\ = 0 and the other one to \Q\/M = 0.2. 


and r/ > —4a. For a < 3 


subject to the constraints r^ > r, 

have a naked singularity because r^ < r m j n . When 3 
there is a critical value of charge now given by 


— 1 + y / 1 — 16MA/(37 t) /(8A), we 
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The number of horizons are two when \Q\ < \Q C \, one (degenerate) for \Q\ = \Q C \ and zero if 
\Q\ > \Q C \- When a > M/tt we have again a naked singularity. As it was previously said the value 
of ao is taken greater than r m i n and r^. 

The critical value of charge for fixed a is an increasing function of A, then for A < 0 its value is 
smaller than for A = 0. For small values of |A|, which are the most interesting ones from a physical 
point of view, a numerical calculation shows that the stability regions slightly change compared 
with the case A = 0. For the sake of brevity, the plots are not included. 


3.3 De Sitter case 


The case A > 0 is considerably different from the point of view of the character of the horizons in 
the original manifold: for A > 0 a cosmological horizon exists. Therefore the shell should always 
be placed inside the cosmological horizon. The positions of the singular surface and the horizons 
are obtained again from Eqs. El and El respectively. When 0 < A < 37r/(16M) we have that 


> T}i for 3 |^—1 — \J\ — 16MA/(37t) /(8A) < a < 3 —1 + y / 1 — 16MA/ (37r) /(8A), so there 


is a naked singularity. With other combinations of the parameters it is r m i n < r/i, if r/,. exists. For 
M/tt — 1/A < a < M/tt, besides the cosmological horizon, there are two, one or zero additional 
horizons for \Q\ < |Q c |i IQ I = IQc| and \Q\ > \Q C \, respectively, with the critical value of charge 
again given by Eq. El - When a < M/tt — 1/A or a > M/tt there are no horizons and we have 
again a naked singularity. As in the other cases, the value of ao is taken greater than r m i n and 77 1 . 

As in the preceding case, |Q C | increases with A, so that the critical value in De Sitter case is 
greater than for A = 0. Now there is a change in the values of a for which the largest regions 
of stability are found, namely M/tt — 1/A < a < M/tt. However, the stability analysis for small 
A does not reveal any remarkable aspect differing from the case A = 0, except that the stability 
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region corresponding to (3 2 < 0 becomes limited for large values of the wormhole radius due to the 
presence of the cosmological horizon. Again, for the same reason as in the Anti-De Sitter case, the 
plots are omitted. 


4 Energy conditions and exotic matter 

Quantifying the amount of exotic matter has been considered as a way to characterize the via¬ 
bility of traversable wormholes [7]. Here we shall analyze the energy conditions and evaluate the 
total amount of exotic matter for the wormholes constructed in Section 2, in the case of static 
configurations, i.e. for a = ao- 

The weak energy condition (WEC) states that for any timelike vector it must be T^u^u 1 ' > 0; 
the WEC also implies, by continuity, the null energy condition (NEC), i.e. that for any null vector 
k^ it must be T^ lu k fJ 'k 1 ' > 0 J2]. In an orthonormal basis the WEC reads p > 0, p + pj > 0 Vj, 
while the NEC takes the form p + pj >0 Vj. In the case of the wormhole constructed above with 
radial pressure p r = 0, we have cr < 0, a + p r < 0, so that both energy conditions are violated. 
The transverse pressure is pt = p and the sign of a + pt, instead, is not fixed, but depends on the 
values of the parameters. 

There have been several proposals for quantifying the amount of exotic matter in the wormhole; 
two of them are the integrals um 

J pV^gd 4 x, j (p + pi)^gd 4 x, (26) 

where g is the determinant of the metric tensor. The most usual choice is the integral including 
the pressure associated to the violation of the energy conditions: 

Q = J (p + p r )V-~9 d 4 x. (27) 


In our case, introducing the new coordinate 77 = ±(r — ao) with ± corresponding to each side of 
the shell, we have 

2n 7T 7T oo 
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III! (p + Pr)\/—9 dlZdO d'x dtp. 


(28) 


0 0 0 —oo 


Because the shell does not exert radial pressure, and the energy density is located on a surface, so 
that p = <5(77) ao, then we simply have 


27T 7T 


n = 


a 


V^d\ r=ao dddxdip = 47r 2 ago- 0 . 


ooo 


(29) 


Thus we find that 

Q = Sira^y/ f(a 0 ). (30) 

Replacing the explicit form of the metric simple expressions for the behavior of 17 with the 
wormhole radius can be obtained for the limiting case a —> oo, which makes sense only for A < 0 
and Aa > —4/3. For A = 0 we have 

I7«-37rag. (31) 
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Figure 6 : Amount of exotic matter (for A = 0) as a function of the charge, for different values of 
the parameter a and the wormhole radius a o- Left: ao = 2.1 right: ao = l.lrv In both cases, 
a\ = O.lM/n, a .2 = O.hM/ir, and 0:3 = 0.9M/7T. The horizon radius depends on the value of a. 
Note the difference in the scales in the energy axis between the two graphics. 


For A < 0 the limiting expressions differ depending on the sign of cc; for a > 0 we obtain 
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while for a < 0 we have 




37raQ 



(33) 


A natural question is which are the conditions such that the amount of exotic matter can be 
reduced. Because hi is proportional to ao, which, as anticipated in Section 2, approaches to zero 
when the wormhole radius tends to the event horizon of the original metric, we shall analyze this 
limit in detail. We consider values of a for which horizons effectively exist; the number of horizons 
in the original manifold is in this case determined by the value of the charge. When \Q\ is less that 
the critical value of charge \Q C \ defined in Section 3, there are two horizons (within this analysis 
we are not interested in the cosmological horizon which appears in the case A > 0 , which is always 
much greater than the others for small A). In this case /( 77 J = 0 and f'{rh ) 7 ^ 0. For ao near the 
event horizon radius rh, we have 
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+ O(a 0 - r h ) 1/2 , 


(35) 


so that as ao approaches rh the total amount of exotic matter 17 tends to vanish, but the pressure 
takes unlimitedly large values. Instead, when \Q\ = \Q C \, f{rh) = ) = 0 and f"(rh) / 0, we 

have 
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4^/2 nr h 


Vf''( r h)(a 0 - r h ) + O(a 0 - r h ) 2 , 
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so that the amount of exotic matter Q can be made as small as desired, keeping the pressure finite, 
by taking ao near rh■ One can also verify that if \Q\ > \Q C \ then the original metric includes no 
horizons, and the amount of exotic matter has a minimum for a q slightly greater than r m i n , while 
the pressure for this ao remains finite. 

In Figs. El we have plotted the total amount of exotic matter for different choices of the param¬ 
eters (always with a > 0 and A = 0), in the case \Q\ < \Q C \, so that at least one horizon exists. 
In each figure, when a changes we keep fixed the relation between the wormhole radius and the 
horizon; the range of possible values of charge is reduced as a increases. We can see that: 1) for a 
given value of a and a fixed radius of the wormhole, the total amount of exotic matter decreases 
as the charge is made larger; 2) for a given charge, the amount of exotic matter is reduced by 
increasing the value of a; 3) the exotic matter present is reduced by placing the wormhole throat 
nearer the horizon. 

5 Discussion 

We have studied five dimensional spherically symmetric thin-shell wormholes in Einstein-Maxwell 
theory with the addition of a Gauss-Bonnet term. We have analyzed the mechanical stability 
under perturbations preserving the symmetry, and evaluated the total amount of exotic matter 
and related its behavior with that of the pressure of the shell, as a function of the parameters of 
the model and of the theory. 

For null a our paper extends to five dimensions the analysis made in previous works HUES] 
in four dimensional gravity. But when a / 0 new interesting results are obtained. We have found 
that the inclusion of the quadratic correction (Gauss-Bonnet term) allows for much more freedom 
in the choice of the configurations to render them stable. In particular, a central positive feature of 
the inclusion of the quadratic correction is that, differing from the general relativity case (a = 0), 
considerably larger regions of stability appear, even for vanishing charge; in this sense, we should 
emphasize that values of the parameter /3 2 positive and smaller than unity are now possible, while 
in the case a = 0 this could be achieved only with the aid of large values of charge (Q close to 
V3M/n). 

Regarding the amount of exotic matter Q and pressure, the analysis shows that, in the case that 
an event horizon exists in the original manifold from which we started our thin-shell construction, 
can be minimized by choosing a wormhole radius near the event horizon. In general this is 
correlated with an unlimited increase of the pressure, except if the parameters are chosen so that 
the inner and outer horizons coincide: in this case the exotic energy can be reduced as desired, 
keeping the pressure finite. Besides, we have shown that as the Gauss-Bonnet term is made more 
relevant by increasing the parameter a, the amount of exotic matter results to be substantially 
reduced. 
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